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The simulation of quantum processes is a key goal for the 
grand programme aiming at grounding quantum technologies 
as the way to explore complex phenomena that are inacces- 
sible through standard, classical calculators ^ [3] . Some in- 
teresting steps have been performed in this direction: simple 
condensed matter and chemical processes have been imple- 
mented on controllable quantum simulators Jsi 01] . The rel- 
ativistic motion and scattering of a particle in the presence 
of a linear potential has been demonstrated in a trapped-ion 
quantum simulator Q-lzl] that opens up the possibility to the 
study of quantum field theories Isl [3- The quantum Ising 
model fid] has been experimentally analysed under the per- 
spective of universal digital quantum simulators and very 
recently scaled up to hundreds of particles |[l2|] . This scenario 
has recently been extended to open quantum evolutions lO . 
marking the possibility to investigate important features of the 
way a quantum system interacts with its environment, includ- 
ing the so-called sudden death of entanglement induced by a 
memoryless environment 

This interaction can destroy the most genuine quantum 
properties of the system, or involve exchange of coherence be- 
tween system and environment, giving rise to mernory effects 
and thus making the dynamics "non-Markovian" lllSn . Char- 
acterizing non-Markovian evolutions is currently at the cen- 
tre of extensive theoretical and experimental efforts 
Here we demonstrate experimentally the (non-)Markovianity 
of a process where system and environment are coupled 
through a simulated transverse Ising model (see Ap- 
pendix). By engineering the evolution in a fully controlled 
photonic quantum simulator, we assess and demonstrate the 
role that system-environment correlations have in the emer- 
gence of memory effects. 

The paradigmatic description of quantum open dynamics 
involves a physical system S evolving freely according to the 
Hamiltonian and embedded in an environment E (whose 
free dynamics is ruled by H^). S yste m and environment in- 
teract via the Hamiltonian H^^ lilSn . which we assume to 
be time-independent for easiness of description. While the 
dynamics of the joint state psE is closed and governed by 
the unitary operator Ut = e~^^^ +^ +^ ^\the state of S, 
which is typically the only object to be accessible directly, is 
given by the reduced density matrix ps{t) = T^^E{psE{t)}- 
For factorized initial states p5'£;(0) = p5'(0)(8)p£;(0) and mov- 
ing to an interaction frame defined by the free Hamiltonian 
of the total system, such reduced evolution can be recast in 
the operator- sum picture ps{t) = (0)^/1' where 



{K^} is the set of trace-preseving, non-unitary Kraus opera- 
tors of S that are responsible for effects such as the loss of 
populations and coherence from the state of the system 

The roots of such decoherence mechanism have long been 
studied, together with their intimate connection to the so- 
called measurement problem and the implications of the col- 
lapse of the wave function 1 19I-I21]. Numerous experimentally 
oriented techniques have been proposed to counteract deco- 
herence Yet, a complementary viewpoint can 
be taken, where the possibility to engineering structured en- 
vironments and tailored S-E couplings is seen as a resource 
to achieve longer coherence times of the system (2?], prepare 
entangled states, perform quantum computation, and realise 
quantum memories. This calls for the exploitation of memory- 
effects typical of a non-Markovian dynamics as a useful tool 
for the processing of a quantum state. 

Unfortunately, a satisfactory understanding of non- 
Markovianity is yet to be reached, which motivates the recent 
and intense efforts performed towards the rigorous character- 
isation of non-Markovian evolutions, the formulation of cri- 
teria for the emergence of non-Markovian features, and the 
proposition of factual measures for the quantification of the 
degree of non-Markovianity of a process lEsl - EsI] . Some 
of them have been recently used in order to characterize, 
both theoretically and experimentally, the character of system- 
environment interactions for few- and many-body quantum 
systems 

Here we use a photonic setup to simulate a system- 
environment coupling ruled by a transverse Ising spin model 
(see Appendix). Such interaction gives rise to the non- 
Markovian evolution of S, as witnessed by the non-monotonic 
behavior of its entanglement with an ancilla A that is shielded 
from the environmental effects. Our goal is the experimental 
investigation of the fundamental connection between system- 
environment correlations and non-Markovianity. The simula- 
tor that we propose allows for the implementation of various 
free evolutions of S and E, as well as the adjustment of their 
mutual coupling, thus making possible the transition from 
deeply non-Markovian dynamics, all the way down to a fully 
forgetful regime. It is realised using different degrees of free- 
dom of the photonic device shown in Fig.[TJ which consists of 
the concatenation of two Mach-Zehnder interferometers and 
a single Sagnac loop. The information carriers are two pho- 
tons (referred hereafter as "high" and "low"): the system S 
is embodied by the polarisation of the low photon (which can 
be either horizontal \H) = |0) or vertical \ V) = |1)), while 



Q HWP @ 22,5° Q HWP @ 0° 

QST ^ 1^) modes , 




FIG. 1 : The Sagnac interferometer used to perform our simula- 
tion. The red line represents the high photon, whose polarisation is 
immediately detected. The low photon, depicted by the yellow and 
blue paths, goes through the set of gates acting on its polarisation 
and momentum, implemented by BS, HWP, and glasses plates 0^. 
The optical axis of HWP* is kept free so as to implement the last 
step of the simulation. The thin glass plates ^i, 02 and 03 placed 
in the interferometer allow to set the phase of the environment evo- 
lution at each step. The Hadamard gates are implemented by setting 
01 = 02 = 03 = [cf. Eq. (□]. 
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FIG. 2: Scheme of principle of the simulation, a) A maximally 
entangled state of A and S is prepared. S and E then evolve through 
the application of free single-qubit gates (H^ and H^) and joint 
two-qubits (d-tkTrCZkTv) ones. Time increases from bottom to top 
of the figure. Here 7Vh{i) stands for the polarization of the high (low) 
photon, while ki for the momentum of the low photon, b) Each block 
of the scheme is technically realised by the concatenation of a BS and 
two HWPs oriented at different angles. 



the environment E is encoded in the longitudinal momentum 
degree of freedom (the path) of the same photon (which will 
be right |r) = |0) or left |/) = |0)). The ancilla A is embodied 
by the polarisation of the high photon. High and low photons 
are emitted by the source of polarization-entangled described 
in Appendix. 

The ancilla is the key tool for our goals. Indeed, to investi- 
gate the emergence of non-Markovianity in the evolution of S 
due to its interaction with E, we use the method proposed in 
Ref . 1I26I1 : we focus on the modifications induced by the S-E 
coupling on a prepared entangled state of S and A. If the S-A 
entanglement decays monotonically in time, the dynamics of 
S is fully Markovian. Differently, if for certain time- windows 
there is a kick-back from E that makes such entanglement in- 
crease, the dynamics is necessarily non-Markovian. In fact, 
if the local action of the environment is no longer represented 
by a continuous family of completely positive maps, the S-A 
entanglement is no longer constrained to decrease monotoni- 
cally. This is evidence of the flow-back of coherence on the 
system and results in an increment of the S-A entanglement. 

As in other digital quantum simulators, the dynamics is 
approximated by a stroboscopic sequence of quantum gates. 
Conceptually, the simulation consists of the forward evolu- 
tion of S over discrete time slices 1I29I1 according to a Trotter- 
Suzuki decomposition fsO^ of the total time propagator (see 
Appendix). This approach is known to be effective for quan- 
tum simulation 0] and is implemented here by the sequence 
of operations shown in Fig. [2^). The free evolution of the 
environment is accounted for by the Hadamard gate = 
{cr^ -\- cr^)/ {criri is the m = x^y^z Pauli matrix of qubit 
j = S^E and W is the Hadamard gate of qubit j = S^E), 
while (X af . The S-E interaction is engineered by im- 
plementing the controlled-rotation G^^ = \0) {0\-^ ^ ^ 
|1) 1Z^{(p), which rotates the system according to the 
general single-qubit operation {(p) = cos (pt^ — i sin (pay 



depending on the state of the environment. This class of 
conditional operations is obtained from Hamiltonian gener- 
ators of the two-qubit transverse Ising form (see Appendix), 
which motivates our choice and is thus the class of system- 
environment interactions that is simulated in this work. Exper- 
imentally, we fix = 7r/4 so as to implement a conditional- 
(a^l-L^) gate. By using the identity af = H^a^H^, it is 
straightforward to prove that this gate is locally equivalent (via 
erf ) to the composition of a controlled Hadamard and a con- 
trolled anti-Z (i.e. a gate applying af only when E is in |0)), 
namely = Cnk^CZk^ = |0) (O|^0af + |1) 
The extra af needed to make the two conditional gates equiv- 
alent can thus be absorbed in , which thus becomes . 

This set of gates is experimentally realised in our photonic 
simulator as sketched in Fig lit). The gate is implemented 
by means of a beam splitter (BS) in conditions of temporal and 
spatial indistinguishability of the optical modes. This scheme 
allows to evolve the input modes |r) and \£) at each step as 
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(Ir)-e^^M^)), (1) 



where the phases 0^ are varied by rotating thin glass plates 
intercepting one of the optical modes entering the BS. In order 
to perform a more general rotation of E it will be necessary 
to unbalance the output modes (using intensity attenuators) 
making the probability of occurrence of |r) and \£) unequal. 

The controlled gate CHk-K is realised by placing a half- 
wave plate (HWP) on one of the two output modes with op- 
tical axis at 22.5° with respect to the vertical direction. The 
temporal delay introduced by this waveplate is compensated 
by another HWP on the opposite output mode with optical 
axis set at 0° with respect to the vertical direction. This im- 
plements the CZ gate, as shown in Fig.[2j3). In order to as- 
sess non-Markovianity, ancilla and system are prepared in the 
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maximally entangled state = {\HH) + \VV))/V2. 

This two-qubit Bell state is engineered using the polarization 
entanglement source of Ref. 113 ill (see Appendix). The en- 
vironment is initialised in |x)f = (a |r) + Vl — oP' 
(a G M is set by the transmittivity of [cf. Fig.[Tl and can 
be varied by placing an intensity attenuator on one of the two 
output modes), a state endowed with quantum coherence as it 
is key due to the conditional nature of the dynamics that we 
simulate. The stability, modular structure and long coherence 
time of our interferometric setting allows for the repeated it- 
eration of each block of gates {U^ ®t^). 

As A does not evolve with the environment, the polarisation 
state of the high photon (depicted in red in Fig.[T]) is immedi- 
ately detected and only the low photon (in yellow in Fig. [T]) 
goes through each sequence of gates at the various steps. A 
standard optical setup for the performance of quantum state 
tomography (QST) [ 32] is used to reconstruct the state of the 
S-A system after the evolution. The radiation is collected 
by using an integrated system composed of GRaded INdex 
(GRIN) lens and single-mode fibre 11331] . and is then detected 
by single-photon counters. For each step of the simulated dy- 
namics, we measure the state of the environment by project- 
ing the polarisation qubits on the states \HH), \HV), \VH), 
\VV). Finally, we trace out the degrees of freedom of the S- 
A system by summing up the corresponding counts measured 
for every single projection needed for the implementation of 
single-qubit QST. The Pauli operators for E, generated after 
the first passage through the ^^i, are measured using BS2- 
The same procedure is followed to perform the QST of the 
environmental state at each step. The state of the S-A system 
is reconstructed in a similar way, by summing up the counts 
collected after projecting E onto |r) and \£). 

Fig. [3] summarises the results obtained by running through 
the evolution of the overall system. In order to quantify entan- 
glement we use the entanglement of formation EOF(SA) ll34l] 
between S and A, which is operatively linked to the cost of en- 
gineering a given state by means of Bell-state resources (see 
Appendix). We are also interested in the correlations shared 
by the environment with the rest of the system, hence we eval- 
uate the von Neumann entropy of the environment, defined as 
S{E) = —Tt[pe log2(p£:)], which under the assumption of 
pure total ASE state quantifies the entanglement in the parti- 
tion AS against E. Figure [3^) shows the experimental entan- 
glement of formation at each time step (black- square points) 
against the results of a theoretical model (red line) that, in- 
cluding all the most relevant sources of imperfections, devi- 
ates from the ideal picture sketched above. First, the BSs are 
not entirely polarisation-insensitive: for BSi and H (V) polar- 
isation the reflectivity over transmittivity ratio R/T = 42/bS 
(45/55), while for BS2 is R/T = 45/55 (55/45). Second, 
although the desired input state 10+)^"^ is created with high 
fidelity (~ 93%, see Appendix), the entangled- state source 
generates spurious \HV) and \ VH) components, accounting 
for about 5% of the total state, which reduce the initial system- 
environment entanglement to about 0.8. The inclusion of such 



imperfections makes the agreement between theory and ex- 
perimental data very good up to the fourth step of our simu- 
lation, showing at least one revival of the S-A entanglement 
and thus witnessing the non-Markovian nature of the evolu- 
tion 1E6I1 . The fifth experimental point is significantly far 
from the theoretical behaviour because of the not-ideal set- 
ting of the phases 0^ (i = 1,2,3). In fact, we have verified 
that their values are not completely polarization independent. 
This determines a slight difference for the four contributions 
\HH)g^, \HV)s^,\VV)s^,\VH)s^ entering the state. This 
imperfection affects the performance of each Hadamard gate 
and becomes significant expecially for the last step, where 
the cumulative effect of three gates should be considered. 
Nonetheless, the last point reveals successfully the occurrence 
of a second entanglement revival, thus strengthening our con- 
clusions. 

Figure [3t)) compares the experimentally inferred EOF(SA) 
(black squares) to the von Neumann entropy of the environ- 
ment Se (green circles) quantifying the correlations shared 
between E and S-A. These figures of merit appear to be 
perfectly anti-correlated, thus giving evidence of a trade- 
off between the amount of entanglement that S and A can 
share at the expenses of S-E correlations. This strength- 
ens the idea that correlations with the environment play a 
fundamental role in this process, a point that has been ad- 
dressed in Ref. |[35ll where it is shown that the establishment 
of system-environment correlations is a necessary conditions 
for the emergence of non-Markovianity. This result can be 
bridged with our analysis considering that the mixed system- 
environment state of Ref. [35] can be purified by enlarging the 
Hilbert space of the system including an appropriate ancilla. 

To reinforce this point even further, we now explore the 
Markovian counterpart of our simulation by replacing the uni- 
tary evolution of E with an incoherent map that resets the en- 
vironment into the very same state at each step of the evo- 
lution. As the key role in the S-E interaction is played by 
quantum coherence, we chose to re-set the environment into a 
completely mixed state. This is realised by spoiling the tem- 
poral indistinguishability of the optical modes entering the 
BS by mutually delaying the |r) and \i) components. Intu- 
itively, by making the state of the environment rigid, we wash 
out any possibility for system-environment correlations, thus 
pushing the dynamics towards Markovianity. This intuition 
is fully confirmed by the experimental evidences collected 
via QST: the black squares in Figure [3]:) show a monotonic 
(quasi-exponential) decay of EOF(S A) (matching our theoret- 
ical predictions, red line), which is in perfect agreement with 
the absence of S-E quantum correlations as signalled by the 
positivity of the partially transposed S-E state. 

By simulating a non-trivial two-qubit coupling model, we 
have demonstrated the non-Markovianity of the evolution 
induced on 6* by a dynamical environment and a system- 
environment interaction allowing for kick-back of coherence. 
By adopting a witness that makes use of the effects that non- 
Markovianity has on entanglement, we have explored experi- 
mentally the link between the emergence of non-Markovianity 
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FIG. 3: Results of the simulation of the (non-)Markovian regimes, a) Theoretical prediction (red line) against experimentally inferred 
data (black squares) for EOF(SA) as against the steps of the non-Markovian simulation. Two revivals of EOF(SA) are clearly visible, b) 
Comparison between the experimental evolution of EOF(SA) (black squares) and S(E) (green circles) representing between E and S-A. The 
two evolutions are clearly anti-correlated. In order to measure the Pauli operators needed to reconstruct the state of E at the fifth step, a further 
BS will be necessary, c) Theoretical behaviour (red line) and experimental data (black squares) for EOF(SA) in the Markovian simulation. 



and system-environment correlations. The next step in this en- 
deavour will be the experimental proof that non-Markovianity 
can be used as a resource for the advantageous process- 
ing of information, such as the preparation of interesting 
states, along the lines of previous studies on state engineer- 
ing and information manipulation through Markovian pro- 
cesses H^IItI]. Our setup will be particularly well suited for 
this task, in light of the effective control over both system and 
environment that can be engineered both in space (thanks to 
the modular nature of the setting) and in time. 
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APPENDIX 



path entangled, as explained in Ref . Bill . For our pourposes 
we need to select a pair of correlated directions, belonging to 
the emission cone surface, along which two photons travel, 
obtaining in this way the high and low photons shown in red 
and yellow in Fig. 1. 



Entanglement of formation 



For two qubits, the entanglement of formation is given 

by 134] EOF = h (^ ^+^^-^' ^, with = -xlog^x-il- 

x) log2(l — x) and C = max{0, Ai — A2 — A3 — A4}. Here, 
{Ai} is the set of eigenvalues (arranged in non-increasing or- 
der) of the Hermitian matrix R = \f~^Jpp^ where p = 
{py (jy)p*{(jy (g) (J y) p* is thc conjugatc of the density 
matrix. 



Source of entangled states 



ControUed-rotation gates from quantum Ising models 



The black box shown in Fig. 1 of the main text represents 
the source used for the generation of the hyperentangled 
state encoded in the polarization and path degrees of freedom 
of the photons produced by spontaneous parametric down- 
conversion. Here polarization entanglement, corresponding 
to the state = -^{\HH) + \VV)), was realized by 

spatial and temporal superposition of photon pair emissions 
occurring with equal probability, back and forth, from a 
single barium borate (BBO) type-I crystal under double 
excitation of a vertically polarized UV CW pump beam. 
A suitable rotation of the polarization on one of the two 
possible emission directions of the photons was then applied. 
The source that we have used is such that the two photons, 
belonging to the degenerate BBO emission cone, are also 



Here we sketch the procedure for the achievement of a con- 
trolled rotation from a quantum Ising model and determine 
the conditions to achieve a desired rotation angle. For conve- 
nience of notation, we will use the correspondences (already 
introduced in the main text) {|0) , \l)}s a ^ {\H) , \V)}s a 
and{|0),|l)}^^{|r),K)}^. _ 

The transverse Ising-like model lliol] whose associated 
propagator has been simulated in our work is 



jjl sing 



■ E 

p=x,y,z 



(2) 



We take ^ J, so that the dynamics of the E subsystem as 
induced by this coupling is effectively frozen. In the compu- 
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tational basis, this model can be written as 



+ \1){1\e® 



(3) 



The diagonal terms proportional to ef only give rise to a shift 
in the energy of the two subspaces that have been identified by 
this formal splitting. We can thus safely neglect them. On the 
other hand, each non-diagonal matrix acting in the S subspace 
gives rise to a non-trivial single-qubit evolution of the form 



cos(z^o^) 



^(ef + J) 



sin(i/o^) 



sin(i/o^) 



— sin(z/o^) 



M) + %±^sinM). 



cos(i/it) 



^ sin(^it) 
— sin(z/it) 



cos(i^it) + - 



■ sin(i^it) 



■ sin(i^it) 
(4) 



where = 2^^(ef )2 + [ef + + (ef )2 (k = 

0,1). Needless to say, we want Uq = l*^ for a time r 
which the conditional gate is performed. On the other hand, 
we should enforce Ui = Tl^{(p) at r. The first condition is 
straightforwardly fulfilled for z/o'?"/(27r) = n G Z. As for the 
second condition, using the definition of IZ^ {(f) given in the 
main text, we should enforce 



cos(p = cos(i^ir), sincp = (e^/i^i) sin(i/iT), 
= [(ef - J)/^i] sin(zyir) = -(ef /z^i) sin(z/ir). 



(5) 



These can be solved, at a chosen value of J, for ef = J, which 
+ 4 and z^i = 27re^. In turn, these give 



give uo 
us 



27r 



27rr' 



47r2 



n. 



(6) 



As the value of J is determined by the experimental setting 
chosen to implement the gate, these conditions fully specify 
the values of the parameters in the Ising model needed in order 
to implement a controlled rotation gate with an arbitrary value 
of the angle cp. 



Experimental system-ancilla density matrices 

In Fig. [4] we report the S-A states that have been created 
experimentally at each step of the non-Markovian dynamics 
that we have simulated. We have used the operative approach 
put forward in Ref. Js^l, which is the standard in the tomo- 
graphic reconstruction of multi-qubit states, to determine the 
density matrices of the S-A system, which have then been ex- 
pressed in the Bell basis spanned by = -i=(|00) ± |11)) 

and|V)* = ^(|01)±|10)). 



[10 
[11 

[12; 

[13 
[14 
[15 

[16; 

[17 

[18 
[19 

[2o; 

[21 

[22; 

[23 

[24; 

[25 
[26; 



Electronic address: " a. chiuri @ uniroma 1 . it I 

Electronic address: l.mazzola@qub.ac.uk 

Buluta, L, & Nori, F. Quantum simulators. Science 326, 108 111 

(2009). 

Cirac, J. I., & Zoller, P. Goals and opportunities in quantum 
simulation. Nature Phys. 8, 264 (2012). 

Ma, X.-S. , Dakic, B., Naylor, W., Zeilinger, A., & Walther, 
P. Quantum simulation of the wavefunction to probe frustrated 
Heisenberg spin systems. Nature Phys. 1 399 (2011). 
Lanyon, B. P., ah Towards quantum chemistry on a quantum 
computer. Nature Chem. 2, 106 (2010). 

Gerritsma, R., et al. Quantum simulation of the Dirac equation. 
Nature 463, 68 (2010). 

Gerritsma, R., et al Quantum Simulation of the Klein Paradox 
with Trapped Ions. Phys. Rev. Lett. 106, 060503 (2011). 
Blatt, R., & Roos, C. F. Quantum simulations with trapped ions. 
Nature Phys. 8 277 (2012). 

Casanova, J. et al. Quantum Simulation of Quantum Field The- 
ories in Trapped Ions. Phys. Rev. Lett. 107, 260501 (2011). 
Casanova, J., Mezzacapo, A., Lamata, L., & Solano, E. 
Quantum Simulation of Interacting Fermion Lattice Models in 
Trapped Ions. Phys. Rev Lett. 108, 190502 (2012). 
Sachdev, S. Quantum Phase Transitions (Cambridge Univ. 
Press, 2011). 

Lanyon, B. P., al. Universal Digital Quantum Simulation with 
Trapped Ions. Science 334, 6052 (2011). 
Britton, J. W., et al. Engineered two-dimensional Ising inter- 
actions in a trapped-ion quantum simulator with hundreds of 
spins. Nature 484, 489 (2012). 

Barreiro, J. T. et al. An open- system quantum simulator with 
trapped ions. Nature 470, 486 (2011). 

Almeida, M. P., et ah. Environment-Induced Sudden Death of 
Entanglement. Science 316, 579 (2007). 

Breuer, H-P. & Petruccione, F. The Theory of Open Quantum 
Systems (Oxford Univ Press, 2007). 

Liu, B.-H., et al. Experimental control of the transition from 
Markovian to non-Markovian dynamics of open quantum sys- 
tems Nature Phys. 7, 931 (2011). 

Apollaro, T. J. G., et al. Memory-keeping effects and forgetful- 
ness in the dynamics of a qubit coupled to a spin chain. Phys. 
i?^v. A83, 032103 (2011). 

Nielsen, M. A., & Chuang, I. L. Quantum Computation and 
Quantum Information (Cambridge Univ. Press, 2000). 
Zurek, W. H. Decoherence and the transition from quantum to 
classical. Phys. Today 44, 36 (1991). 

Zurek, W. H. Decoherence, einselection, and the quantum ori- 
gins of the classical. Rev Mod. Phys. 75, 715 (2003). 
Schlosshauer, M. Decoherence, the measurement problem, and 
interpretations of quantum mechanics. Rev. Mod. Phys. 76, 
1267 (2005). 

Zanardi, P., & Rasetti, M. Noiseless Quantum Codes. Phys. Rev. 
Lett. 79, 3306 (1997). 

Viola, L., & Lloyd, S. Dynamical suppression of decoherence 
in two-state quantum systems. Phys. Rev. A 58, 2733 (1998). 
Myatt, C. J., et al. Decoherence of quantum superpositions 
through coupling to engineered reservoirs. Nature 403, 269 
(2000). 

Wolf, M. M., Eisert, J., Cubitt, T. S., & Cirac, J. I. Assess- 
ing Non-Markovian Quantum Dynamics. Phys. Rev. Lett. 101, 
150402 (2008). 

Rivas, A., Huelga, S. F. & Plenio, M. B. Entanglement and 



6 



Non-Markovian Dynamics 




FIG. 4: Tomographic reconstruction of the state shared by S and A at each step of the Non-Markovian dynamics. The initial state is represented 
by a 1 0^ )^^. It is worth noting that in the fourth step the coherence terms are null, in fact in this case we obtained the lowest value of EOFsa- 
In the fifth step there is a revival of the coherence terms but with smaller values with respect to the third step. The imaginary matrices are 
negligible for each reconstructed state. 



non-Markovianity of quantum evolutions. Phys. Rev. Lett. 105, 
050403 (2010). 

[27] Breuer, H-R, Laine, E-M., & Piilo, J. Measure for the degree 
of non-Markovian behavior of quantum processes in open sys- 
tems. Phys. Rev. Lett. 103, 210401 (2009). 

[28] Lu, X.-M., Wang, X., & Sun, C. P. Quantum Fisher Information 
Flow in Non-Markovian Processes of Open Systems. Phys. Rev. 
A 82, 042103 (2010). 

[29] Lloyd, S. Universal Quantum Simulators, Science 273, 1073 
(1996). 

[30] Trotter, H. F. On the product of semi-groups of operators. Proc. 

Am. Math. Soc. 10, 545 (1959). 
[31] Barbieri, M. et al. Polarization-momentum hyperentangled 

states: Realization and characterization. Phys. Rev. A 72, 

052110(2005). 



[32] James, D., Kwiat, P., Munro, W. & White, A. Measurement of 
qubits. Phys. Rev A 64, 052312 (2001). 

[33] Rossi, A., Vallone, G., Chiuri, A., De Martini, F, & Mataloni, 
P. Multipath entanglement of two photons. Phys. Rev. Lett. 102, 
153902 (2009). 

[34] Wootters, W. K., Phys. Rev Lett. 80, 2245 (1998). 

[35] Mazzola, L. Rodriguez-Rosario, C., Modi, K., & Paternostro, 
M. Dynamical role of system-environment correlations in non- 
Markovian dynamics, arXiv: 1203.3723 (2012). 

[36] Kraus, B. et al. Preparation of Entangled States by Quantum 
Markov Processes. Phys. Rev A 78, 042307 (2008). 

[37] Vertraete, R, Wolf, M. M., & Cirac, J. I. Quantum computa- 
tion, quantum state engineering, and quantum phase transitions 
driven by dissipation. Nature Phys. 5, 633 (2009). 



